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1. INTRODUCTION 


The computer program used to determine the rigid and elastic stability 
derivatives presented in the summary reoort (Ref. 1) is listed in this appendix 
along with instructions for its use, sample input data and answers. This program 
represents the airplane at subsonic and supersonic speeds as (a) thin surface(s) 
(without dihedral) composed of discrete panels of constant pressure according to 
the method of Woodward (Ref. 2) for the aerodynamic effects and slender b^am(s) 
for the structural effects. 

Given a set of input data, the computer program calculates an aerodynamic 
influence coefficient matrix [A] and a structural influence coefficient mat-ix [C] 
by the methods specified above. In a different option [CJ can be read in from tape. 

The stability derivatives determined by this program and the equations used 
in their calculation are given in Table 1. From the table it can be seen that for 
rigid stability derivatives (as well as AC p ) (A] is needed, whereas for elastic 
stability derivatives both [A] and [C] are required. 

The program is compatible with the geometry requirements and definitions 
of reference 3, using assumptions defined In the Appendix to this report. 

Using the method of Reference 4 the program is also capable of computinq 
C D*/ C L Z f° r rigid and elastic wings. 
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2. TABLE OF SYMBOLS 


The units used for the physical quantities defined in this paper are given both in 
the International System of Units (SI) and the U.S. Customary Units. 

Input Data 


Symbols 

Dcscr i pt ion D Imens ion 

ALT 

Altitude 

Feet (m.) 

ALTNUM 

Number of altitude 
runs to be made 


AM 

Mach Number 


AMASS (I) 

Mass of panel “1" 

slugs (Kgm) 

[c] 

Structural influence 
coefficient matrix 

rad/lb. (rad/N.) 

c.. 

Structural influence coefficient. 

rad/lb (rad/N) 

•J 

angle of attack induced on panel 
i due to a unit load on panel j 


CLDES 

Desired lift coefficient 


CPCWL 

Constant percent chord- 
wise lines 


CPSWL 

Constant percent stream- 
wise lines 


CREF 

Reference Chord 

Feet (m.) 

CSHT 

Structural root chord of 
the horizontal tail 

Feet (m .) 

CSMW 

Structural root chord of 
the main wing 

Feet (m.) 

El (I) 

El value of the Ith elastic 

lb - inch^ (N-m^) 


axis segment 

or lb - ft^ 

GJ (I) 

GJ value of the 1th elastic 

lb - inch? n 


axis segment 

or lb - fr (N-rn) 

IASIGN (K) 

The Kth element of the 
array IASIGN and the number 
of an elastic axis end-point to 
which the Kth panel is attached 


KON1L 

= 1 for fuselage geometry data 
= 2 for wing geometry data 



= 3 for horizontal tail geometry data 
= 4 for canard geometry data 
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Symbols 

Description 

Dimension 

LPNL (1) 

Number of panels on front 
fuselage 


LPNL (2) 

Number of panels on rear 
fuselage 


LPNL (3) 

Numbers of panels on main 
wing 


LPNL (4) 

Number of panels on horizontal 
tail 


M 

Maximum number of elastic 
axis end-points 


MOVTEL 

= 0, Non-movoble toil 
= 1, Movable tail 


N 

Totol numbei of panels or, 
total number of unit loading points 


N1 

Attachment point of horizontal 
tail elastic axis to rear fuselage 
elastic axis, point number 


NC 

Number of sections of the 
aerodynamic surface 


NCNTL 

See Sec. 4 for complete description 


NCONT 

= 0 if geometry data are input 
according to the present program 

= 1 if geometry data are prepared 
accordinq to Ref. 1 


NCONTL 

= 1 for front fuselage data 
= 2 for rear fuselage data 
= 3 for main wing data 
= 4 for horizontal t. 1 -lata 


NCP1 

Number of CPCWL on fuselage 


NCP2 

Number of CPCWL on wing 


NCP3 

Number of CPCWL on horizontal tail 


NCP4 

Number of CPCWL on canard 


NCPSWL 

Number of constant pe -ent stream- 
wise lines 


NEA 

Number of elastic axis end-points 
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Dimension 


Symbol 5 

NE1 GJ 

RHO 

SOS 

SREF 

TITLE 

W 

WDTHHT 

WDTHMW 

XCG (1) 
XEm (I) 
XREF 

XXL (1) 

XXL (2) 

XXT (1) 

XXT (2) 

YCG (1) 

YEA (1) 

YL (1) 

YL (2) 

Z 


Description 


= 0 if El and GJ values are in 

lb-inch^ 

= 1 if El and GJ values are in 
ib-ft 2 


Air density 

Slugs/ft'* (Kg/m 2 ) 

Speed of sound 

ft/ sec . (m/ sec) 

Reference area 

Any title for computer run 

ft 2 (m 2 ) 

Total weight of the airplane 

lb. (N.) 

One-half the width of the fuselage 
at the horizontal tail 

Feet (m) 

One-half the width of the fuselage 
at the main wing 

Feet (m) 

x-coordinate of the loading 
point "I" 

Feet (m) 

x-coordinate of the elastic 
axis end-point "I" 

Feet (m) 

Attachment point of the main 
wing elastic axis to the fuselage 
elastic axis, x-coordinate 

Feet (m) 

x-coordinate of the L.E. of 
inboard chord 

Feet (m) 

x-coordinate of L.E. of 
outboard chord 

Feet (m) 

x -coordinate of T .E . of 
inboard chord 

Feet (m) 

x-coordinate of T.E. of 
outboard chord 

Feet (m) 

y-coordinate of the loading 
point 1 

Feet (m) 

y-coordinate of the elastic 
axis end-point "I" 

Feet (m) 

y-coordin Jte of inboard chord 

Feet (m) 

y-coordinate of outboard chord 

Feet (m) 

Z-coordinate of the section under 
consideration 

Feet (m) 
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Output Data 


Symbols 

Description 

Dimension 


Aerodynamic influence coefficient 
matrix 

ft^rad ^ (m^rad 

a.. 

'J 

Aerodynamic influence coefficient, 
load on panel i resulting from a unit 
angle of attack on panel j 

ft^rad ^ (nwad 

ALTITUDE 

Altitude 

Feet (m) 

AREA 

Area of a panel 

(Feet) ^ (m^) 

CD 

Total induced drag coefficient 


CDI 

Sectional induced drag coefficient (Cj.) 


CL 

Total lift coefficient 


CLALPA 

, Variation of lift coefticient 
a g with angle of attack for the 

elastic case with zero mass 

per rad. 

CLALPE 

C L , Variation of lift coefficient 

flp with angle of attack for the 

elastic case including mass 
effect 

per rad. 

CLALPR 

C, / Airplane lift curve slope 
L a 

per rad . 

CL1 

Sectional lift coefficient (C|) 


CLQ 

, Variation of lift coefficient 
q with pitch rate 

per rad. 

CLQEBR 

C. , Variation of lift coefficient 

with pitch rate for the elastic 
case with zero mass 

per rad. 

CLQI 

C L / Inertial ly induced variation of 
qj lift coefficient with pitch rate 


CLTDDI 

C, , Inertial ly induced variation 

of lift coefficient with pitch 
angular acceleration 

2 

sec. 

CLTRIM 

c, 

Lt _* 
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Symbols 

CLWD1 

CMALPA 

CMALPE 

CMALPR 

CMQ 

CMQEBR 

CMQ1 

CMTDD1 

CMWDI 

CP 

CT 

DCLDN 


Description 


Dimension 






Inertial ly induced variation 
of lift coefficient with rate 
of downward velocity per- 
turbation 

Variation of pitching 
moment coefficient with angle 
of attack for the elastic case 
with zero mass 

Variation of pitching 
moment coefficient with angle 
of attack for the elastic case 
including mass effect 

Variation of pitching moment 
coefficient with angle of attack 
(i.e., static longitudinal stability) 

Variation of pitching moment 
coefficient with pitch rate 

Variation of pitching moment 
coefficient with pitch rate 
for the elastic case with 
zero mas 


C , Inertial ly induced variation of 
^q. pitching moment coefficient with 

pitch rate 

C Inertial ly induced variation of 

pitching moment coefficient with 
pitch angular acceleration 

Z , Inertial ly induced variation of 
m w. pitching moment coefficient with 

rate of downward velocity 
perturbation 

AC p , Change in pressure coefficient 

between upper and lower surfaces 

Sectional leading edge thrust coefficient 


sec. per ft. 

(sec^/ m ) 

per rad. 
per rad . 

per rad. 

per rad . 
per rad. 


2 

sec. 

2 r*. 

sec . per ft 

(sec.V m ) 


dC. /dn , Variation of lift coefficient with 
load factor 
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Symbols 

Description 

Dimension 

DCMDN 

dC /dn , Variation of pitching 
moment coefficient 
with load factor 


DYNAMIC PRESSURE 

Dynamic pressure 

lb per ft 2 (N/m 2 ) 

Eld) 

El Value for the Ith elastic 
axis segment 

lb— ft 2 (N-m 2 ) 

GAMMA 

Non-dimensional local 
circulation (Cj c /2b) 


GJ (I) 

GJ value for the 1th elastic 
axis segment 

lb- ft 2 (N-m 2 ) 

IASIGN (I) 

t urrber of elastic axis end- 
point to which panel "I" is 
attached 


MASS (I) 

MaSs of Ith panel 

Slugs (Kg) 

REF. AREA 

Reference area 

(Feet) 2 (m 2 ) 

REF. CHORD 

Reference chord 

Feet (m) 

RKO 

Air density 

Slugs/ft 2 (Kg/m 2 ) 

SOS 

Speed of sound 

ft. per. sec. (m/sec) 

VELOCITY 

= (Mach Number) x 
(speed of sound) 

ft. per. sec. (m/sec) 

WEIGHT 

Weight of the airplane 

lb. (N) 

XCG (I) 

v- coordinate of the centroid 
of Ith panel or the Ith unit 
loading point 

Feet (m) 

XCP (1) 

x-coordinate of the control 
point of 1th panel 

Feet (m) 

YEA (1) 

x-coordinate of the 1th 
elastic axis end-point 

Feet (m) 

Y 

Non-dimensional spanwise 
location 


YCG (I) 

y-coordinate of the centroid 
of Ith panel or the Ith unit 
loading point 

Feet ( m ) 

YCP (l) 

y-coordinate of the control point 
»i Ith panel 

Feet (m) 

YEA (I) 

-coordinate of the Ith elastic 
uxis end-point 

Feet (m) 

ZCP (I) 

z-coordinate of the control pomt.^_ . 

of Ith panel 

Feet (m) 
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3. PREPARATION OF INPUT DATA 


The preparation of input data for this computer program is done in three steps 
Step I Preparation of planform geometry data (Sec. 3.1) 

Step 2 Preparation of mass distribution data (Sec. 3.2) 

Step 3 Preparation of structural data (Sec. 3.3) 
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3.1 Preparation of Planform Geometry Data 


The ground rules for defining the planform geometry and for dividing a planform 
into aerodynamic panels are described below. 

A general planform is taken as an example. Its projections in the x-y plane and 
the x-z plane are given in Figure 1 . 

The follov/ing paneling definitions should be observed. 

1 Break Lines 

A break line on a planform is a stream line which connects the leading and trailing 
edges and occurs when either the leading - or trailing - edge has a slope discontinuity. 
Root and tip chords are also considered as break lines. A break line on one lifting surface 
creates a break line on the others. (See Figure 1 .) 

Referring to the x-y plane view of the planform, there are 14 break lines. Table 2 
defines these break lines precisely. 


TABLE 2 


Break Line 


Explanation 


1- r 

2 - 2 ' 

3- 3' 

4- 4' 

5- 5' 

6 - 6 ' 

7- 7* 

8 - 8 ' 

9- 9' 

10 - 10 ' 
11-11’ 
12 - 12 ' 

13- 13' 

14- 14* 


Rout chord of fuselage 

Tip chord of fuselage 

Root chord of wing 

Due to fuselage-wing overlapping 

Due to break line 12-12' on horizontal tail 

Due to slope discontinuity at 6' 

Due to break line 14-14' on horizontal tail 

Due to slope discontinuity at 8 

Tip chord of wing 

Root chord of horizontal tail 

Due to horizontal tail-fuselage overlapping 

Due to slope discontinuity at 12 

Due to break line 6-6' 

Tip chord of horizon*-! tail 


2 Sections 

A section is a region on the configuration between two consecutive break lines. For 
the example planform of Figure 1 there are seen to be 11 sections. 

The variable NC on input card 4 sets the number of sections. 

3 Constant Percent Chordwise Lines (CPCWL) 

These are the lines which divide ihe sections in the chordwise direction. In general 
there could be any number of CPCWL's in o section. In the current program it is limited to 
35. These can be different or different aerodynamic surfaces. The 0% and 10^% lines a r e 
the leading and trailing edges of the section. For accurate values of Cq / , these 

lines are selected by using the scheme given in Ref. 3. i 

4 Constant Percent Streamwise Lines (CPSWL) 

These are the lines which divide the section in the streamwise direction. In general 
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there could be any number of CPSWL's In a section. In the current program It Is limited 
to 35. However, in case one section is behind another section in the sfreomwise direction, 
or one section Is directly above another section, it Is desirable to use the same CPSWL's. 
The 0% and 100% lines are the inboard and outboard chords of the section. 

5 Panel Numbers 

Panels are numbered in the Increasing order of section numbers, from inboard to out- 
board of the section and from the leading edge to the trailing edge. 

6 Total Number of Panels (N) 

The present computer program is limited to 300 panels. 


7 Control Point of a Panel 

In all cases the downwash control point is located at 0.95 of the local chord which 
passes through the centroid of the panel . 


13 



3,2 Preparation of Mass Distribution Data 


In the wing-body- tail program there are three different types of mass distributions 
to be considered: wing mass, body mass and tail mass. A procedure for determining 
these masses cannot be given such that all types of airplanes are covered by it. For the 
family of parametric wings studied herein, a detailed procedure for the calculation of 
mass distribution is included in Appendix D of the Summary Report (Ref. 5). It is possible 
to develop such procedures for tails and fuselages of parametric families of airplanes. 

The user of this complete WBH program has the following options: 

a) read in a complete mass matrix from an external source. 

b) use the procedure developed for the parametric family of wings in 
Ref. 1 to find the wing part of the mass matrix. Read in the 
other mass contributions from an external source. 

c) neglect the effect of mass and mass distribution. 

Whether or not option c) is realistic depends on the type of airplane under study. In Ref. 1 
it was shown that in many instances the effect of mass is negligible. For transport type 
(low load-factor type) airplanes this is definitely not so. 

3.3 Preparation of Structural Data 


Given an El distribution for fuselage and El and GJ distribution for wing and hori- 
zontal tail, the method given In Appendix E of the Summary Report (Ref. 6), is used to 
determine all the input data for structural matrix. 
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This computer program has five different options, depending upon the value of NCNTL , which 
appears on input data card number 1 . 
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Number of sections of the aerodynamic surfaces. 
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If CCREF ~ 0 Program takes the Inboard chord of the last 
processed region as reference chord. 
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ALT Altitude (ft.) 

Cord 22 is repeated "ALT MUM" times. 



5. OUTPUT DATA FORMAT 


The output data are divided into five sections. 

1 . Geometry data 

The planform geometry data is outputted exactly in the same way as it was 
inputted. The format is self-explanatory in the output printout. The output vari- 
ables are: 

a) The number of sections into which the planform is divided. 

b) For each section inboard and outboard chord coordinates are given along with 
the constant percent chordwise lines and constant percent streamwise lines for 
that section. 

c) REF. AREA and REF. CHORD represent the reference area and the reference 
chord respectively. 

2. Rigid longitudinal derivatives 


Co mputer Variable Explanation 

CLALPR r. 

L a 

CMALPR c 

m « 

CLQ C L 

9 

CMQ r 


3. Sectional C^/Cy^ and the total induced drag parameter Qy/CL^ for rigid 
configuration 

Computer Variable Explanation 


Y 

Non-dimensional spanwise location 

CDI 

Sectional induced drag coefficient (C,j;) 

CLI 

Sectional lift coefficient (C^) 

GAMMA 

= CjC/2b 

where C - local chord 

CT 

Sectional leading edge thrust 

CDI/CL* * 2 

* C di /C|2 

CD/(CL ' cl; 

Total Induced Drag Coefficient 

C, 2 
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4. ACp values for ora= l .0 radians 

Computer Variable 
PANEL NUMBER 

(XCP, YCP, ZCP) 

(XCG, YCP = YCG) 

AREA 

CP 

MASS 

5. Panel assignment output 

Computer Variable 

I 

(XCG, YCG) 

IASIGN 

6. Elostic axis data 
Computer Variable 

I 

(XEA, YEA) 

El 

GJ 

7. Elostic longitudionol derivatives 

Computer Variable 
WEIGHT 

RHO 

SOS 

ALTITUDE 

VELOCITY 


Explanation 
Panel number 

Control point location 

Panel centroid location 

Panel area 

A Cp 

Mass of the panel 
Explanation 

Panel number or unit loading point 
number 

Panel centroid location or unit loading 
point location 

Number of elastic axis end-point to which 
panel "I" is attached. 

Explanation 

EHstic axis end-point number 
Elastic axis end-point location 
El value of the Ith elastic axis segment. 
GJ value of the Ith elastic axis segment. 

Explanation 

Total weight of the airplane 

Air density and speed of sound at that 
altitude 

= (Mach Number X Speed of sound) 
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DYNAMIC PRESSURE 


Dynamic pressure 


CLTRIM 

c L 

frim 

CLALPA 

c L 

S 

CMALPA 

Cfn r 

CLQEBR 

C L 

q E 

CMQEBR 

c m 

q E 

CLQI 

C L 

q l 

CMQI 

c 

q I 

CLWDI 


CMWD1 

c m . 

w i 

CLTDDI 

c, .. 
*1 

CMTDDI 

c .. 

DCLDN 

dc L/an 

DCMDN 

^ C m/dn 

CLALPE 


CMALPE 

S 
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CDI/C l * * 2 
CD/(CL * CL) 


Explanation 
Panel Number 



Daramete*- Cn:/C| 2 for elastic 


Explanation 

Non-dimensional spanwise location 
Sectional induced drag coefficient (C^;) 
Sectional lift coefficient (C ^ ) 

= qC/2b 
“ C di ./C L 2 


Total Induced Drag Coefficient 



25 



6. TEST CASE 1 


This test case deals wl 
listing Is given in Table 3. 


th Wing 5 of ^ef. ^ (see Fig. 2). An inpur data cards 
Output data for this test case are presented in Table 4. 
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Figure 2. Test Cose 1 Plonform 


Tobte 3. Listing of Input Doto Cords for Test 1 Case 1 
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14.9027 7.59 16,0427 9.40 17.3727 11,43 18.7627 13.67 



Toble 3 (Continued) Listing of Input 1 Dota Cord; for Test Cose 1 


_____ _ ^ 

" REPRODUCIBILITY OF THE ORIGINAL PaGE IS POOR. 
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Table 4. Output for Test Case 1 

TtST CASS 1 
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7. TEST CASE 2 


This test case deals with a complete elastic wing, body and horizontal tail configur- 
ation. Figure 3 defines the overall geometry for this case. Figure 4 defines the elastic 
axis location. Input data cards are listed in Table 5. Output data listing for this case : s 
given in Table 6. 


38 













2.083FT 
(.6341 M) 







Toble 5. Listing of Input Data Cords for Test Cosa 2 
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Toble 5 (Continued) Listing of Inpui’ Dgfo Cords for Test Cose 2 
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Toble 5 (Continued) Listing of Input Do to Cords for Test Cose 2 
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Table 6, Output for Test Case 2 
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8. TEST CASE 3 


In this test case the geometry data of one of the examples given in Ref. 3 is used. 
This geometry data is reduced by the program to calculated rigid stability derivatives and 


Figure 5 gives the actual geometrical description and geometry assumed 
Input data cards are listed in Table 7 and output data are listed in 
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Figure 5 (Continued) Actual and Assumed Plonform Geometry for Test Cose 3 




Figure 5 (Continued) Actual and Assumed Flapform Geometry for Test Cose ** 
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THIS C *88 IS COD'kSTS 



9. COMPUTER PROGRAM LISTING 


The computer program was originally written by members of the Flight Research 
Laboratory at the University of Kansas for the Honeywell 635 computer system. The 
version listed in this document is a revised version as optimized by Computer Sciences 
Corporation of Hampton, Virginia for the CDC 6600 computer at NASA Langley Research 
Center. 
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FIGURE 6. (continued) 
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FIGURE 6. (continued) 
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FIGURE 6. (continued) 
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READ (5,58) (GARB, I = 1 » NC ANOR ) DRGD2070 

GO TO 26 DRGD2080 

25 NCANOi =-NCAN0R DRGD2U9G 

READ ( 5 , 58 ) ( GARB * I = 1 ♦ NCANOI ) DRGD2100 



RFAD (5»58)(GARB*I=1» NC AN01 ) DRGD2110 

RFAD ( 5,58) (GARB, I=1»NCAN01 > DRGD2120 

26 CONTIMJF - DRGl>2 130 
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IF (NCAN.FO.l) GO TO 33 DRGD2420 
XFUM3)=XHT(?»?) DRGD2A30 
YRIM3 )=YHT(2,2 ) DRGD2440 
ZFUM3 )=C HRDHT(2,? ) DRGDZ450 



ZAVHT ( n = (ZHT( 1,1 )+ZHT( 1 ,2 >+ZHT <2 >2)1/3. DRGD2A60 


oooo oooooo ooooooooooooooooooo 

f^coOO-^ • r\l cO 3“ lT\ sO i v '00 0'OHMm<riN'ONa)(^0^(Njf04‘ir) 

u> l ninxo>o>o>OsOO^OONONOr^i s ^r'r^r^r^ 
(Mrsjrgf\jrNjrsjc\jfvif\jrg <MrMfNrM(\j'sj(\icgr\jr\jr>jrNj'MrMr\/r\jr\jcsjrvj 
OOQQOOQOQO OQQOQOOCiOQCi^OCiOOQQa 
OOOOOi/vOOOO ooooooooooooooooooo 
csraaacraoraraa: ccoc&cnocx.acx'. xiccxctcxiccctccQccza' 
0030000000 OQOOOOOQOOOOOOOOOOO 
I I 

cu cu 


tO 


Z It 
< < 



3 


3 


UJ 








U- 


ll 


u 








> 


>- 


< 








■w 


w 


Ll 








V, 


N 


3 X 








— 


— » 


2 3 








X 


X 


</> 








h~ 




1 - 












Ll O 








in 


m 










3 


o 


-J *-» 








u. 


u_ 


>— 








1 


1 


UJ LL 










— * 


2 — 










r — 1 


O -1 








H.H 


- 










l-t 


i — 1 


2 Q£ 




— * 




in 


m 


O UJ 




CM 


— 


3 


3 


X 




2 


1—4 


iL 


Li- 


m h- 




2 


Ll 


V 


> 


LU o 


— * 


*■ 


< 


w> 




2 


r-W 


#-< 


3 


* 


* 


*-* |jj 


2 


2 


LU 


— 


— * 


3 X 


2 


2 


Z 


— 


— - 



r\! 

r 

#-4 

» 


* 


fH 


f-i 


X 


no 


in 


r-f 




— 


< 2 

-l 

•> 

3 

3 

r*t 

z 


f -4 


*— 1 


UJ o 

M 

«-H 

H 4 

U- 


z 


^n 


•n 


c: 

< 

in 

< 

X 

3 

*> 




3 

• 

cc m 

V— 

3 

H- 

» 

HI 

m 


Ul 


Ll 

CM 

LU 


Ll 


CM 

< 

3 


X 


NJ 

X 

H 2 

O 

>* 

o 

w- 


Ul 


1 


1 

*— * 

3 ~ 

k- 

•> 

h- 

rH 


X 






O 3 


CM 


in 

c 

* 


J 


(M 

CM 


UJ 


LU 

3 

H- 

Ul 





» 

m v 

3 

in 

3 

Ll 


< 


#— < 


f— * 

(VJ 

o < 

O 

3 

O 

rs! 

UJ 

* 


m 


in 

«— 

2 UJ 


IL 


•> 

3 

UJ 


3 


X) 

K 

M Q£ 

r -1 

rsl 

ro 

*— • 

C 

z 


Ul 


Ll 

X 

Ll CD 




lA 


r 

— 

X 


Ni 

tNj 


•J 

lO 

3 

3 

in 

O 

M 

w 



+ 

UJ 


3 

Hi 

Ll 

z 


#> 

1 


1 


X X 

< 

Ll 

< 

X 

•— 

o 

^ — CM 



— 


<t h- 

f- 

>- 

*- 

•> 


X 

, w 

H 


*~4 

» 

UJ 


#• 


r-* 


SJ 

A « 

HH 



X cm 

o c o 

2 

in 

2 

m 

z 


cv cm X 

r— 


H 

i- w 

CD H- 

O 

3 

O 

3 

o 

X 

— — o </) 


in 

3 : »- 



LL 


LL 


#1 

hhQTD 


3 

in X 

UJ 

m 

X 

in 

X 

in 

X 

fMlXIli. 


Ll 

3 M 

UJ 3 

UL' 

•— 

UJ 

H- 

a 

w 

•* X > U X 


Kl 

Ll — 

2 O 

z 

Vf 

z 

vr 

z 

X 

f- II It tl 

II 

H", 

it 

— tl it 


HI 

< 

Hi 

< 


< 

H 



•— 

«pi"% 

H* UJ 

-U 

LU 

3 

LU 

3 

UJ 

*— « H-l • HI 

t— • 

*—* 


i— » (V 

3 U 


CC 


a: 


a: 

— w 

W 

«— 

W 

V 

O < 


cr 

V 

cr 

V 

CD 

r— i |M r— » r— * 

1—1 


r-* 

»— • *— i I— 

or u. 

< 


< 


< 


IT. LT IT, Kf) 

IT m 

m in i 

cc cc 

LU 

3 

UJ 

3 

u 

3 

2 D 3 D 3 3 

33 > 

3 3 

a 

3 

a: 

3 

or 

3 

C U- iL LL 

Ul Ul Ul 

iLlL < 

m m 

X 

< 

cc 

<r 

<r 

< 

CX^NX 

> M 

>- >■ M 



U 


i 


U 


C 

fr 


to 


UUUUOUU u u i*> o o uuu 




92 


DRG0276J 

BRFAK LINES ON WING DUE TO TAIL 2 DRGD2770 

DRGD2700 

CALL BREAK. ( X»Y , CHORD* NEW AF 1 *YFUS1 ,NN2 »NtWAF2 ) DRGD279U 

0RGD2 600 



BREAK LINES ON TAIL 1 DUE TO WING DRGD2810 

DRGD2820 
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HOWEVER. IF AN END OF FILE IS ENCOUNTERED AFTER THE FIRST READ. LN2LU296 

AN ATTEMPT WILL BE MADE TO READ THE ENTIRE MATRIX FROM TAPE7. LN2L0297 
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IF ( NCNTL I 80.80.90 LN2L0327 

C LN2L0328 

C**MATRIX WAS READ BY COLUMNS. STORE BY ROWS. LN2L0329 

80 CALL CWRITE (PHUP( 1 ) .C.NCP. IND) LN2L033U 



90 REWIND 7 LN2L0331 

RETURN LN2L0332 

LN2LU333 
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CALL KFAr>l*S(20*WB»LPANEL. IAA ) LN3L0670 

HO 29 J = 1 »LPANFL LN3...0680 

P(I»JI=WB(J1 N3L069U 

DO 3 1 = 1 vL PANEL LN 3LO lOO 
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PRESSURE COEFFICIENTS (ENDS/ LN3L1U20 

I.N3L1J3-J 

WRITE I 6 *?7 ) LN3L1O40 

WRITE ( 6 *?8 ) (I»WA(I>»CP(I)»1 = 1 »L PANEL ) LN3l_loS0 



CALL SET DRAG (XN*YN*YCP*AREA»CP»XCG»WA»CL»P»N) LN3L1U60 

18 CONTINUE LN3L107U 

RETURN LN3i_10flO 
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10. Appendix 


Assumptions for Data Reduction from Reference 3 


The present computer program is valid only for thin elastic aeroploncs and therefore 

requires the reduction of the configuration data of Reference 3 to satisfy the input data 

format. The following assumptions are made for data reduction . 

1 . The “Y" locations of the chords in Ref. 3 are assumed to be break 
lines for the present program. 

2. The height of the wing (ZAVWNG) is calculated by taking the average value of 
z- coordinates of all break lines. In other words, camber is taken to be zero. 

3. If the fuselage is not circular, it is modified to a circular siape by assuming a 
radius equal to the maximum y- coordinate at each station, the center of which 
is located at the corresponding z- coordinate. 

4. The z- coordinate of the fuselage (ZAVFUS) is found by taking average value of 
centers of circular section. 

5. Half of the fuselage width (FUSWTH) is calculated by taking the average radius of 
fuselage sections at the root chord. 

6. If J(ZAVFUS-ZAVWNG)] < (FUSWTH/2), the wing is assumed to be in the plane 
of the fuselage, i.e., ZAVWNG = ZAVFUS and FUSWTH is replaced by the y- 
coordinate of the root chord of the wing. 

7. If the wing and fuselage are not in the same plane, the wing is extended to the 
center line of the fuselage and the number of break lines for wing is increased by one. 

8. The portion of the fuselage ahead of the wing leading edge is replaced by a 
trapezoid having a width equal to FUSWTH and area equal to the projected area 
of the fuselage in the x-y plane between the nose and wing root leading edge 
(See Figure A1 ). 

9. Assumption 8 is also made for the portion of fuselage behind wing root trailing 
edge. (See Figure Al). 

10. The y- coordinates of a horizontal tail and a canard are assumed to be the y- 
coordi nates of the break lines for the corresponding surfaces. 

11 . The heights of a horizontal tail fzAVHT(l)) and a canard (ZAVHT(2)J are calculated 
by taking the average value of the break line z- coordinates of the corresponding 
surfaces . 

12. The semi-width of the fuselage at the root chord of the horizontal tail (WDTH1) 
and the canard (WDTIT2) is calculated by taking the average radius of fuselage 
sections between the corresponding root chords. 
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13. If I (ZAVHT (I) - ZAVFUS)I < (WDTH 1/2), the toil is assumed to be in the plane 
of the fuselage i.e. ZAVHT(l) = ZAVFS. 

14. If the tail is not in the plane of the fuselage, the tail is extended to the center 
line of the fuselage and the number of break lines of the tail are increased by one. 

15. Assumptions 13 and 14 are also applicable to a canara surface. 

16. A break line on the wing induces a break line at the same y- location on the 
horizontal tail and on the canard and vice versa. 

17. Between the two break lines on any aerodynamic surface the constant percent 
streamwise lines are calculated such that they are at least (b/20) apart. 
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